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We introduce a number of random matrix models describing the Google matrix G of directed 
networks. The properties of their spectra and eigenstates are analyzed by numerical matrix diago- 
nalization. We show that for certain models it is possible to have an algebraic decay of PageRank 
vector with the exponent similar to real directed networks. At the same time the spectrum has no 
spectral gap and a broad distribution of eigenvalues in the complex plain. The eigenstates of G are 
characterized by the Anderson transition from localized to delocalized states and a mobility edge 
curve in the complex plane of eigenvalues. 

PACS numbers: 89.75.Hc, 89.20.Hh, 89.75.Fb 


I. INTRODUCTION 

The phenomenon of Anderson localization [I] appears 
in a variety of quantum physical systems including elec¬ 
tron transport in disordered solids and waves in random 
media (see e.g. It is usually analyzed in the frame 

of Hermitian or unitary matrices. Recently, the localiza¬ 
tion properties of nonunitary complex matrices has been 
analyzed for Euclidean matrices [4] in relation to light 
and wave localization [5]. 

In this work we analyze the possibilities of Anderson 
like localization and delocalization for the matrices be¬ 
longing to the class of Markov chains and Google matrix 
G [6j [7]. Such matrices have real nonnegative elements 
with the sum of elements in each column being equal to 
unity. For a directed network one first defines an adja¬ 
cency matrix Aij which has element 1 if a node j have a 
link pointing to node i and zero otherwise. The columns 
with only zero elements ( dangling nodes ) are replaced 
by columns with 1/N where N is the matrix size. The 
elements of other columns are renormalized in such a 
way that their sum becomes equal to unity (JA Sij = 1, 
Sij = / JA Aij). Thus we obtain the matrix Sij of 

Markov transitions. Then the Google matrix G of the 
network takes the form El El: 

Gij = QiSij + (1 — ol)/N . (1) 

Here, the damping factor a is taken in the range 0 < a < 
1. In the context of the World Wide Web (WWW) the 
term (1 — a) describes for a random surfer a probability 
to jump on any node of the network. The above con¬ 
struction of G has been proposed by Brin and Page [6] to 
describe the structure of the World Wide Web (WWW). 
For the WWW it is assumed that the Google search en¬ 
gine uses a « 0.85 [7]. We can also consider a generalized 
case of weighted Markov transitions Sij corresponding to 
real positive elements of Aij like happens for the world 
trade network (see e.g. i)- 

The matrix G belongs to the class of Perron-Frobenius 
operators, its largest eigenvalue is A = 1 and other eigen¬ 
values have |A| < a mm- The right eigenvector at A = 1, 


which is called the PageRank (GP = P), has real non¬ 
negative elements P(i) and gives a probability P{i) to 
find a random surfer at site i. It is possible to rank 
all nodes in a decreasing order of PageRank probability 
P(K(i)) so that the PageRank index K{i) counts all N 
nodes i according their ranking, placing the most pop¬ 
ular nodes at the top values K = 1,2,3.... Usually for 
many real directed networks the distributions of number 
of ingoing and outgoing links are described by a power 
law (see e.g. HOD, generating an average approximately 
algebraic decay of PageRank probability P(K) oc 1/K& 
with /? ~ 0.9. Some examples of directed networks can 
be found in mi. 

It is important to note that matrices of Google class 
practically have not been studied in physical systems 
even if they naturally appear in the frame of Ulam net¬ 
works generated by the Ulam method for dynamical maps 
in a coarse-grained phase space (see e.g. EMU). 

Therefore, it is interesting to see if the phenomena of 
Anderson localization and Anderson delocalization tran¬ 
sition can appear for Google matrices. Certain indica¬ 
tions on a possible Anderson transition for the Ulam net¬ 
works, built from dissipative maps, have been reported 
in m with more detailed discussions presented in mi. 
Thus, it would be useful to find random matrix models 
which are able to reproduce typical properties of spec¬ 
trum and PageRank decay in real directed networks. 
However, the results presented in m show that the full 
matrix G with random matrix elements have an unreal¬ 
istic spectrum and hence other random matrix models 
of G should be developed. The models discussed in [IJ| 
give certain indications of delocalization of eigenstates of 
G but the spectrum of G in these models has a large gap 
and is very far from the spectra of real directed networks. 
With this aim we describe below a number of random 
Google matrix models and analyze the properties of their 
spectra and eigenstates. We use certain spectral proper¬ 
ties of small size orthostochastic matrices with N = 3,4 
established in mi- 
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FIG. 1: Google matrix eigenvalues A (a,b), and IPR £ of 
eigenvectors as a function of ReA (c,d). Panels show data for 
RMZ3 model (Sec.2.1) at fixed amplitudes Ei — 0.5 (a,c) and 
for random amplitudes 0.15 < £i < 0.3 (b,d). The green circle 
shows | A | = 1 (a,b); the green horizontal line shows maximal 
possible i — N (c,d); the gray band in (a,c) highlights specific 
states (see text). Here the total number of nodes is iV=8000. 


II. RANDOM MATRIX MODELS OF G 

We start from a description of various random matrix 
models of Google matrix G presenting the results of their 
spectral properties in next Section. 


A. Model RMZ3: random three-diagonal blocks 

Following m we consider ortho stochastic matrix 
blocks Bij of size M x M = 4x4. The orthostochas¬ 
tic property means that = Oij 2 , where an orthogo¬ 
nal matrix O has random matrix elements obtained via 
random rotations. Since O is an orthogonal matrix the 
matrix B is bistochastic with JT Bij = JT Bij = 1 [17] . 
The main reason to use such blocks B is a similarity of 
complex spectrum of random matrix ensemble of B with 
the spectrum of G of university networks of Cambridge 
and Oxford, as discussed in m- The size 4x4 can be 
considered as preferential random links between a group 
of 4 friends. However, a weak point of the random en¬ 
semble of B m is a small matrix size N = 4, while for 
the above universities we have iV^2x 10 5 . 

To go to large values of N in matrix Sij we construct 
the Random Matrix model Z3 (RMZ3) as follows: 
we place blocks B of size M — 4 on the main diagonal 
with weights (1 —ef) and on two adjacent upper and lower 
diagonals with weight Si/ 2 , where Si (i = 1 ,... , N/M) 
are random numbers uniformly distributed in some inter¬ 
val (s m i n , £ max)] each block represents a random realiza¬ 
tion of B ; then the matrix G of total size N is built from 



FIG. 2: Same as in Fig. [l] for RMZ3S model (Sec.2.2) with 
0.15 < Si < 0.3, shortcut amplitude e s = 0.3, <5 = 0.1 (a,c) 
and 5=1 (b,d). Here N = 8000. 


S via the equation 0 . Here we consider two cases with a 
constant Si = 0.5 and the interval range 0.15 < £i < 0.3 
(see Fig. [l]). Obviously, by construction the final matrix 
belongs to the Google matrix class. We use notations Sz 
and Gz for the matrices S and G of this model. 


B. Model RMZ3S: RMZ3 with shortcuts 


The model RMZ3S is obtained from RMZ3 by adding 
shortcut links between blocks B in the upper triangle 
of the whole matrix S , the blocks of shortcut links are 
placed randomly in this part of S. The amplitude of 
transitions from one block to another block (outside of 
three-diagonal blocks of RMZ3) is taken at some fixed 
value e s . The shortcut blocks are randomly and uni¬ 
formly distributed over the upper triangle of the whole 
matrix. After adding the shortcut blocks the columns 
affected by shortcut blocks are renormalized to unity. 
In this way the obtained matrix S again belongs to the 
Google matrix class. The blocks of shortcuts are placed 
randomly in the upper triangle of matrix S, their num¬ 
ber N s is determined by the parameter S = 4 N s / (3 N) . 
In fact 5 gives the ratio of shortcut blocks to the num¬ 
ber of blocks 3N/4 in the model RMZ3. Again each 
block B in the main three-diagonal part of RMZ3 and 
blocks at shortcut positions are taken as random and in¬ 
dependent realizations for each block. We note that the 
shortcuts between single nodes have been used for studies 
of quantum chaos and Anderson transition in the small 
world Anderson model (see miasm The results for 
RMZ3S model are shown in Figs. [5] [3] 
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FIG. 3: Dependence of P(K) for models RMZ3S (Sec.2.2) 
in panels (a,b) and RMZ3F (Sec.2.3) in panels (c,d); here 
0.15 < Si < 0.3. In panel (a) we have 5 = 0.1 (black symbols) 
and S = 1 (blue symbols) at N = 8000; the fitted algebraic de¬ 
pendence is shown by straight dashed lines with parameters: 
a = -6.67, p = 0.288 at 5 m 0.1 and a = -3.76, p = 0.71 
at <5 = 1; panel (b) shows the dependence of P on 5 with the 
full curve for N — 8000 and triangles, crosses and circles for 
N = 2000, 4000,16000 respectively; the amplitude of shortcut 
elements is e s — 0.3. In panel (c) we have 5 — 0.01 (black 
symbols) and 5 = 3 (blue symbols) at /x = 0.1 and N m 8000; 
the fitted algebraic dependence is shown by straight dashed 
lines with parameters: a m —8.39, [3 = 0.072 at 5 = 0.01 and 
a — —6.17, P — 0.36 at 5 = 3; panel (d) shows the depen¬ 
dence of P on 5 for /x = 0.1 (blue) and /x = 0.3 (red) with the 
full curves for N m 8000 and triangles, crosses and circles for 
N = 2000,4000,16000 respectively. 


C. Model RMZ3F: RMZ3 plus triangular matrix 


The results obtained in |T5j 2(1 show that a triangular 
matrix of Google matrix class has a tendency to have a 
realistic PageRank probability decay with P ocl/K and 
have some eigenvalues of finite amplitudes |A|. Due to 
these indications we construct a matrix Sp in the follow¬ 
ing way: N u random numbers fi from the interval (0,1) 
are placed on random positions of the upper triangle of 
matrix of size AT, then all columns are renormalized to 
unity and columns with all zero elements are replaced by 
columns with all elements 1/N. Then we construct the 
matrix G of the model RMZF as: 



FIG. 4: Spectrum (a,b) and IPR £ dependence of ReA for 
the model RMZ3F (Sec.2.3) at 5 = 0.1 (a,c) and 5 = 3 (b,d); 
here 0.15 < £i < 0.3, /x = 0.1, N = 8000; circle and horizontal 
lines are as in Fig. [l] 



FIG. 5: Panel (a) shows dependence of maximal IPR £ (for 
states with |A| < 1) on parameter S for the model RMZ3F 
(Sec.2.3) at N = 8000. Dependence of maximal IPR £ on 
N is shown in panels (b) at 5 = 3 and (c) at 5 = 0.1; error 
bars show statistical error, if it is larger than symbol size, 
obtained from N r disorder realizations. We use N r = 11 at 
N = 2000, N r = 8 at N = 4000, iV r = 4 at N = 8000, 
iV r = 3 at N = 16000. In all panels ji — 0.1 (circles) and 
/x = 0.3 (triangles), 0.15 < Si < 0.3; the straight green lines 
show dependence £ = N; the straight gray lines in (b,c) show 
the fitted dependence (see text). 


D. Anderson models AD2 and AD3 for G matrix 


Szf = (1 —/i) S z ~\~ i^S p , Gzf = qcSzf + (1 — a)/N . (2) 

Here /x determines a measure of contribution of Sp with 
0 < /x < 1. The number of nonzero random elements N u 
is given by parameter 5 = N U /(12N). The results for the 
RMZF model are shown in Figs. |3|4[ [5] 


We use the usual Anderson model mm with diagonal 
disorder terms W{ and transitions V to nearby sites on a 
lattice in dimension d\ 

W\ipi + Vipi+1 + Vipi-1 = Xtpi , (3) 

where indexes in bold are vectors in d-dimensional space. 
On the basis of © we construct the matrices S and G. 
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Thus we consider the dimensions d = 2,3 correspond¬ 
ing to square and cubic lattices. The matrix S is con¬ 
structed as follows: each transition matrix element, cor¬ 
responding to V terms, in the Anderson model in dimen¬ 
sion d Q is replaced by a random number Si uniformly 
distributed in the interval [0, £ ma:c /2d], the diagonal ele¬ 
ment W{ is replaced by unity minus the sum of all Si over 
2 d nearby sites (1 — YalL i e i)- The asymmetric matrix S 
constructed in this way belongs to the Google matrix 
class. Thus we obtain the matrices Sad 2 ,Gad 2 for the 
model AD2 and Sad 3 ->Gad 3 for the model AD3 for 
d = 2 and 3 respectively. The results for these models 
are presented in Figs. 


E. Anderson models AD2S and AD3S with 
shortcuts 


By adding shortcut links between pairs of nodes ran¬ 
domly distributed in the upper triangle of matrix S we 
obtain models AD2S and AD3S respectively from 
models AD2 and AD3. The number of shortcut el¬ 
ements in S is taken to be N s = 2 dN5, their am¬ 
plitude is taken as 0 < Ei < = £max/ 2, after 

adding shortcuts the columns with shortcut elements 
are renormalized to unity. Thus the sum of elements 
in each column is equal to unity and S belongs to the 
Google matrix class. We note the matrices of these 
models as S A D2S , G A D 2 S , Sad3S , Gad 3 S respectively for 
d — 2,3. The results for these models are presented in 


Figs, pm, 10 
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FIG. 6: Distribution of IPR £* on A—plane for the Anderson 
type models AD2 at d = 2 (a) and AD3 at d = 3 (b) (Sec.2.4) 
and the Anderson type models with shortcuts AD2S (c) and 
AD3S (d) at S = 2 (Sec.2.5). Here £m ax m 0.6; N = 130 2 = 
16900 for (a,c); N = 25 3 = 15625 for (b,d) and £max = 0.6, 
e s — Smax /2 = 0.3 for (c,d). Color bars show the ratio £i/N 
(IPR values are averaged inside cells of coarse-grained lattice 
60 x 60). 


we use a = 0.85. The right eigenstates 'ipi(j) of G are 
determined by 


N 


j' = 1 




(4) 


Anderson models AD2Z and AD2ZS with blocks 
and block shortcuts 


By replacing matrix elements in the model AD2 by 
blocks B of size 4x4 (see Sec.2.1) we obtain the model 
AD2Z. In a similar way for the model AD2S we obtain 
the model AD2ZS with block shortcuts. In this case 
we restrict our studies only for dimension d = 2 since 
the matrix size becomes too large for d = 3. Amplitudes 
£max and are defined as for the models AD2 and AD2S. 
Since the transitions are now given by blocks then the 
parameter S is now defined as N s = 2d(N/4)S with d = 2. 
The results for models AD2Z and AD2ZS are presented 
in Figs. Ell m 


We characterize by the Inverse Participation Ratio 
(IPR) & = (Ej ld(j)| 2 ) 2 /Ei |^i(j)| 4 - This quantity 
is broadly used in the studies of Anderson localization 
[3] and determines the number of sites effectively pop¬ 
ulated by an eigenstate. The value of £ is independent 
of normalization. We use normalization Yi P(i) = 1 for 
the PageRank eigenstate at A = 1. For each eigenvector 
we can order all nodes in a monotonically decreas¬ 
ing order of IV^COI thus obtaining the local rank index 
K for a given 'ipi(j). Such a ranking was used in [21J [22] . 
Below we describe the results for models of Sec. 2. 


A. Results for RMZ3 model 


III. SPECTRAL PROPERTIES OF G MATRIX 
MODELS 

We use exact numerical diagonalization for analysis of 
spectrum and eigenstates of models of Sec. 2. The matrix 
size N is changed from a minimal N = 900 up to maximal 
N = 27000. For the description of the decay of PageRank 
probability we use a fit In P = a — (3 In K which gives us 
the PageRank exponent of algebraic p. In all simulations 


For the model RMZ3 at £i = const = 0.5 the spec¬ 
trum is shown in Fig. [lji. We see that it has a form of 
6—rays star typical for the directed networks studied in 

HD EH 122]. The size of the star is slightly reduced since 

all A$(a) a\i(a = 1) for a < 1, except A = 1 0. 
There is also addition reduction of | A* | due to finite cou¬ 
pling terms Si > 0 but this reduction is moderate and the 
spectrum of Gz is close to the spectrum of independent 
4x4 blocks found in H3. Thus RMZ3 model captures a 
part of real properties of directed networks. 
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An interesting property of eigenstates becomes visible 
from the dependence of £ on ReA shown in Fig.[]J at Si = 
0.5. Many eigenstates have relatively small £ < 10 which 
remain bounded with the increase of N up to N = 16000 
(data not shown). However, there is a group of states 
(gray band) with £ ~ N growing linearly with N (data 
not shown). These are delocalized states. Their origin 
becomes clear from the following consideration. We can 
use the anzats in which the elements of 'ip(j) are constant 
inside a given block with a values (p m . Then Eq. @ 
takes the form 

(1 s)(p rn T ^iSPm -\-1 T l)/ 2 — A(^ m , (5) 

since the matrix G is bistochastic with sum of elements 
in rows being unity since Si — const. The spectrum A in 
([5| is real. Thus we obtain in (J5| the Bloch equation with 
plane wave delocalized solutions well known for crystals 
mm- These solutions belong to the gray band part of 
the spectrum in Fig. EK Another part of the spectrum 
corresponds to such 'ip(j) which have different values on 
a scale of one block B. 

The case with different e m (e.g. 0.15 < < 0.3 in 

Fig. & ,d) we can use the same anzats for the left vector 
ip m that leads to the eigenvalue equation: 

(1 ^m)^m T £m( ( tpm +1 T tym— 1)/2 = A(/? m . (6) 

Such a problem corresponds to the case of off-diagonal 
disorder in the Id Anderson model where the localization 
length, and hence IPR, is diverging at the center of the 
band mm- The spectrum A in ( |6|) is real. A similar 
problem is known as the Sinai walk j23j where transition 
probabilities on a Markov chain are fluctuating. This 
model has been studied extensively (see e.g. [24 and 
Refs, therein). 

The spectrum A in corresponding to this anzats, 
is the same for the right eigenvectors [7;. The right 
eigenvectors are different from the left ones but have a 
similar structure on average. The IPR values, shown in 
Fig.p are significantly reduced, comparing to the case 
£ m = const , except those with A close to unity. When N 
is increasing we find that IPR is growing only for A —>> 1 
while for |A| < 1 IPR values remains finite. This corre¬ 
sponds to the known results for the Anderson model with 
off-diagonal disorder. Other eigenstates for which ipi is 
not constant inside B blocks correspond to the eigen¬ 
states with rather small IPR values £ ^ 10. 

Even if the spectrum and eigenstates have interesting 
properties in the two above cases of model RMZ3 there 
is a weak point here: the PageRank probability P in 
these cases is flat being practically independent of K and 
£ ^ N. Thus the situation is very different from the real 
directed networks with (3 ~ 1 (see e.g. [7, H|). This 
happens due to a space homogeneous structure of the 
matrix G (a part of fluctuations) and thus there is no 
leading node with a large number of links. Due to that 
we try to introduce shortcut links as described in the 
next Sec.3.2. 


B. Properties of RMZ3S model 


The spectrum and IPR dependence for RMZ3S model 
with shortcuts are shown in Fig. [2] for two typical values 
of parameter S. We see that at small values of 5 (e.g. 
S = 0.1) the spectrum structure is practically the same 
as for RMZ3 model. However, for larger values (5=1) 
the size of the spectrum star is decreasing. The values 
of IPR are significantly reduced at finite values of 5 and 
our data show that the maximal £ values remain less 
than £ = 200 even for the largest size N = 16000 for 
0.1 <5 < 1 for all |A| < 1 (data not shown). Thus in 
this model all eigenstates remain localized. 

Even if all states are localized the decay of PageRank is 
more close to the case of real directed networks. Indeed, 
the data of Fig. Hf ,b show that P(K ) have approximately 
algebraic decay with PageRank index. The fit allows to 
determine the PageRank exponent /? which is small at 
5 ~ 0.1 and is growing with increase of 5 reaching values 
(3 ~ 0.75 at 5 = 1. It is important to note that (3 is inde¬ 
pendent of N at large N values. Thus the homogeneous 
random elements in the upper triangle of S matrix allow 
to obtain /3 close to unity at large 5. Indeed, in the limit 
of rather large 5 we come to the case of triangular matrix 
S studied in 15] (and also in [2D]) where one obtains an 
approximate decay P oc 1/K. Indeed, at large 5 a sum 
of elements in a row of G drops approximately as 1/K 
(where AMs a row index) leading to P oc 1/K. Indeed, we 
can say that P(K) ~ V • Gxj^j ~ 1/K, where Cj = 1/N 
is a homogeneous initial vector, considering this as one 
iteration of the PageRank algorithm [7]. We note that for 
the PageRank vector we have £ ^ N for (3 < /3 C = 1/4. 

Thus the model RMZ3S has a reasonable spectrum 
structure and an algebraic PageRank probability decay. 
But all eigenstates with |A| < 1 remain localized. Thus 
we go to the analysis of RMZ3F model. 


C. Results for RMZ3F model 

The spectrum and IPR values for the RMZ3F model 
are presented in Fig. [4] We see that the star spectrum 
structure is preserved but IPR values are increased in 
a vicinity ReA « a. The examples of P(K) and (3(5) 
dependencies are shown in Fig. (3j:,d. It is important to 
note that /J is independent of N at large N values. Qual¬ 
itatively, the situation is similar to the model RMZ3S 
but the effect of 5 on localization properties of £ is more 
complicated. 

Indeed, it is well seen in Fig. [5^i that the maximal IPR 
values (excluding PageRank vectors) are at first reduced 
with an increase of 5 from 10 -3 up to 0.1 but they are 
increased when 5 goes from 0.1 to 10. The dependence 
of maximal £ on N at 5 = 0.1; 3 is shown in Fig. Et’C 
for /i = 0.1; 0.3. We fit this dependence by a power law 
£ oc N 1 ' and obtain for /i = 0.1: v = 0.352 (at 5 = 0.1) 
and v = 0.770 (at 5 = 3); for fi = 0.3: v = 0.33 (at 
5 = 0.1 and 3). These results show that there are certain 
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states (except PageRank) that become delocalized in the 
limit of large matrix size. In a certain sense, for the 
dependence £(7V) we have a certain similarity with the 
results obtained in m where a sub-polynomial growth 
of £ with N has been found for randomized university 
networks and preferential attachment models. However, 
for the RZ3F model the spectrum has no large gap and 
is more similar to the real directed networks. 

The investigations of RMZ3F model at larger sizes (e.g. 
with the help of the Arnoldi method JTl] [2T] ) can provide 
a more firm conclusion about the delocalization proper¬ 
ties of eigenstates in this model. 


D. Properties of AD2 and AD3 models 



The spectra of AD2 and AD3 models are shown in 
Fig. §t ,b with color plot of IPR values. We see that 
there are rather large values of £ indicating existence of 
delocalized eigenstates. Indeed, a more detailed analy¬ 
sis presented in Fig. [7] shows that for the states of the 
spectral range ReA > 0.25 IPRs are growing with N 
clearly demonstrating delocalization. Indeed, for max¬ 
imal £ from this range (excluding PageRank) we find 
v = 0.75 at d = 2 and v = 0.95 at d = 3. At the 
same time in a vicinity of the spectrum edge ReA <0.25 
we have v = 0.18; 0.05 for d = 2; 3 clearly showing that 
in this part of the spectrum the eigenstates are well lo¬ 
calized. Indeed, for these localized states we have an ex¬ 
ponential decay In |^| oc —K l / d with the eigenstate rank 
index K (see Fig. [9k, b). Such a decay also appears for 
the localized states of the Anderson model in dimension 
d. 


But for the majority of eigenstates we have significant 
growth of £ with N showing that these states are delo¬ 
calized. Of course, the case of d = 2 should be studied 
in more detail since for the standard Anderson model at 
d = 2 Q all eigenstates are exponentially localized [3]. 
However, we have here non-Hermitian matrix and for our 
knowledge there are no rigorous results about localization 
in such matrices in d = 2. 

Even if in AD2, AD3 models we find delocalization, 
the PageRank in these models is practically flat due to 
absence of central node (see Fig. [§i ,b). Another weak 
point of AD2, AD3 models is a relatively narrow distri¬ 
bution of eigenvalues with |ImA| <0.1 and due to that 
we continue our analysis with the next model. 


E. Results for AD2S and AD3S models 

The spectra of AD2S, AD3S models are shown in 
Fig. ,d. We see that the additional terms in upper 
triangle of matrix S produce a broadening of ImA which 
however still remains relatively narrow (|ImA| < 0.2). 
The IPR values are growing with N except of the eigen¬ 
states at the spectral edge ReA ~ 0.6 (see Fig. [8|. For 
these localized states the exponent v is practically zero 


FIG. 7: Dependence of £ on ReA (a,c) and £ on A (b,d) for 
the models AD2 (a,b) and AD3 (b,d) (see Sec.2.4). For AD2: 
panel (a) is for N = 900 (blue, N r = 10 realizations) and 
N — 16900 (red, N r = 1); panel (b) shows dependence £(A) 
with fits £ oc N 1 ' for eigenstates at the spectrum edge with 
ReA = 0.23 — 0.25 (triangles, v — 0.18) and for maximal £ 
(circles, v — 0.75). For AD3: panel (c) is for N = 1000 (blue, 
N r = 10 realizations) and N = 15625 (red, N r = 1); panel 
(d) shows dependence £(A) for eigenstates at the spectrum 
edge with ReA = 0.23 — 0.25 (triangles, v — 0.05) for maximal 
£ (circles, v — 0.95). Here £max = 0.6. The fits are shown 
by gray lines, green (b,d) and blue, red (c,d) dashed lines 
show dependence £ = N. For panels (b,d) the number of 
realizations changes from N r — 10 to 3 when N changes from 
minimal to maximal value. 


while for the maximal IPR (except PageRank) we find 
rather large values of v = 0.57 at d = 2, v == 0.73 at 
d = 3. Thus, in these models we clearly have the Ander¬ 
son type transition from localized to delocalized eigen¬ 
states. 

In analogy with the 3d Anderson model [3], we make a 
conjecture that in models AD2, AD3, AD2S, AD3S there 
is a certain mobility edge curve in the complex plane A 
which separates localized from delocalized states. In a 
qualitative manner such a curve is visible in Fig. [6] as a 
border between blue color of localized states with small 
£ and white color of states with large £. But definitely 
more detailed studies are required for a more exact de¬ 
termination of such a mobility edge curve. 

Examples of PageRank probability decay are shown 
in Fig. [9] The new element, appearing in AD2S, AD3S 
models (comparing to AD2, AD3 cases), is a dependence 
of the PageRank exponent /? on the parameter S as shown 
in Fig. [lO} These data demonstrate that /3 increases from 
/3 « 0.2 at 5 = 0.1 up to /3 « 0.9 at 5 = 3. Thus 
AD2S,AD3S models have delocalized eigenstates and the 
PageRank exponent of real directed networks. The only 
weak point is a narrow distribution of spectrum in ImA. 
To improve this feature we study in next Section the 
models AD2Z, AD2ZS. 
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FIG. 10: Dependence of the PageRank exponent on the 
parameter 5 for the models AD2S (left panel) and AD3S (right 
panel). Left panel: the solid curve shows data for N — 80 2 , 
with triangles for N = 40 2 and circles for N = 130 2 . Right 
panel: the solid curve shows data for N m 20 3 , with triangles 
for N — 10 3 and circles for N — 25 3 . Here £max = 0.6 and 
= 0.3 


FIG. 8: Same as in Fig. [7] but for the models AD2S (a,b) 
and AD3S (c,d) (see Sec.2.5) at 5 = 2; all parameters are 
as in Fig. m The fits give: (b) v = 0.04 at the spectrum 
edge around Re A w 0.6 (triangles), v — 0.57 for maximal £ 
(circles); (d) v — —0.19 at the spectrum edge around Re A ~ 
0.6 (triangles), v — 0.73 for maximal £ (circles). Here £max — 
0.6, e s — £max/2 — 0.3. For panel (b) [(d)] the number of 
realizations changes from N r = 10 to 3 [1] when N changes 
from 900 to 16900 [27000]. 



FIG. 9: Dependence of eigenvector amplitudes |^| on their 
rank index K for models AD2 (a), AD3 (b) from (Sec.2.4) 
and AD2S (c), AD3S (d) from (Sec.2.5). Here 5 = 0 for (a,b) 
and 5 — 2 for (c,d); N = 10 4 for (a,c) and N — 20 3 for 
(b,d). We use £max = 0.6 and = 0.3 in (c,d). Data show 
maximally delocalized (maximal £ corresponding to PageR¬ 
ank, magenta upper curve) and maximally localized (smallest 
£, blue bottom curve) eigenstates. 


F. Results for AD2Z and AD3ZS models 


The spectra of AS2Z, AD2ZS models are shown in 
Fig. [lT] We see that the star structure appears due to 
introduction of blocks 4x4. The dependence of IPR £ on 
ReA clearly shows the existence of two groups of states 
with small £ < 100, presumably for localized phase, and 



FIG. 11: Spectrum A (a,b) and IPR £ vs. ReA (c,d) for 
the models AD2Z (a,c) and AD2ZS at 5 = 0.25 (b.d) from 
Sec.2.6. Here iV = 4 x 70 2 = 19600, Smax = 0.6 and = 0.3 
in (c,d). 


large £ > 100, presumably for delocalized phase. 

The distribution of £ on A—plane is shown in 
Fig.[l2fc,b. Again we see signs of the mobility edge curve 
separating localized (blue) and delocalized (white) eigen¬ 
states. 

The dependence of £ on A is shown in Fig.[l2fc. There 
are well localized states with £ practically independent of 
N (£ < 20) and delocalized states for which £ is growing 
with N with a relatively large growth exponent v = 0.67 
at 5 = 0 and v = 0.53 at 5 = 0.25. This gives a strong 
argument for existence of the Anderson transition with a 
mobility edge in a complex A—plane in these models. 

The decay of PageRank probability is shown in 
Fig. 12 i: at 5 = 0 we have a flat P(K ) distribution 
with the exponent (3 — 0.16, while at 5 = 0.25 we find 







































FIG. 12: Top panels show distribution of IPR £ values on 
A—plane for models AD2Z (a) and AD2ZS at S = 0.25 (b) of 
Sec.2.6 with parameters of Fig. EH color bar gives the ratio 
£/iV obtained from cells as in Fig. 6] Panel (c): dependence 
of £ on N for AD2Z with triangles for states with A located 
in the delocalized domain ReA £ (0.3, 0.85) (red triangles, fit 
gives v — 0.67) and in the localized domain ReA < —0.5 
(blue triangles, v — 0.15); for AD2ZS at 6 = 0.25 with 
circles for states with A located in the delocalized domain 
ReA £ (0.2,0.85) (red circles, v — 0.53) and in the quasi- 
localized domain ReA < —0.5 (blue circles, v — 0.25); fits are 
shown by lines, green line shows £ = N. Panel (d): depen¬ 
dence of PageRank probability P on PageRank index K for 
models AD2Z (gray symbols) and AD2ZS at S = 0.25 (black 
symbols); the fits for the range K £ (100, 6000) are shown by 
dashed lines with f3 — 0.16 (AD2Z) and f3 — 0.51 (AD2ZS) 
for the parameters of panels (a,b). 


/? = 0.51 being close to the values found in real directed 
networks (e.g. for the Twitter network /3 ~ 0.54 DU)- 
Thus we can say that the model AD2ZS is the one 


being most close to real directed networks with the num¬ 
ber of interesting features: algebraic decay of PageRank 
probability with the exponent (3 ~ 0.5, absence of spec¬ 
tral gap at a = 1, a broad star like distribution of eigen¬ 
values in the complex A—plane, existence of localized and 
delocalized eigenstates of the Google matrix with strong 
indications on the Anderson transition and the mobility 
curve in A—plane. 

We expect that a similar model AD3ZS constructed in 
dimension d = 3 from the AD3S model will have even 
stronger delocalization properties. 


IV. DISCUSSION 

In this work we described various random matrix mod¬ 
els of the Google matrix of directed networks. Our results 
show that for certain models (like AD2ZS) we have an al¬ 
gebraic decay of PageRank probability with the exponent 
P ~ 0.5, absence of spectral gap at a = 1, existence of 
the Anderson transition and mobility edge in the com¬ 
plex A— plane. We think that the further analysis of the 
models described here will allow to establish more close 
links between the Anderson delocalization phenomenon 
in disordered solids and delocalization of eigenstates of 
the Google matrix of directed networks. 
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